Search for: All records

This review focuses on how the modeling of dense granular media has advanced over the last 15 years. The jumping-off point of our review is the μ( I) rheology for dry granular flow, which opened the door to generic flow field modeling but was primarily geared toward problems involving small monodisperse grains of simple shapes. Our review focuses on advances in modeling more material types and behaviors including new approaches for modeling finite-grain-size effects or nonlocality, polydispersity and unmixing, and nontrivial grain shapes. We also discuss growing application areas with tractable order-reduction strategies with a focus on intrusion and locomotion problems. 

We study how channel width variations influence the dynamics of free-surface granular flows. For this purpose, we extend a continuum model framework to granular flows passing through channels that narrow or widen. Our theory uses a linearized approximation to an established dense granular flow rheology and a Coulomb friction law to model interaction between flow and sidewalls. We test the theoretical predictions using two novel 40 cm-diameter drums (convex and concave) filled halfway with 2 mm diameter particles rotated at rates in which the shear layer remains shallow and dense. We apply particle tracking velocimetry to enable quantitative comparisons between experimental data and theoretical predictions. We find that our experimental kinematics and energy profiles largely agree with the theoretical predictions. In general, flows through narrowing channels are faster and deeper than flows through widening channels. The influence of width variations grows with increasing flow speed, and the form of the rate dependence changes fundamentally as the regime changes from one in which kinetic energy is dissipated locally to one in which it is advected downstream. For both regimes, theoretical scaling analysis leads us to experimentally validated power laws, in which the exponent depends on the flow regime, and the multiplicative coefficient depends on channel geometry alone. Finally, we discuss how the differences between theoretical predictions and experimental data may be useful for improving our understanding of flows through non-uniform channels. 

Abstract Alluvial fans are large-scale depositional structures commonly found at the base of mountain ranges. They are relatively soil-rich compared to the rocky terrains, or catchment areas, from which their material originates. When frequented by debris flows (massive, muddy, rocky flows) they contribute significantly to local hazards as they carry focused, collisional, fast-moving materials across alluvial fans, unpredictable in size, speed, and direction. We research how fine particle content in debris flows correlates with directional changes, i.e., debris flow avulsions. Toward this, we analyzed field data from two neighboring alluvial fans in the White Mountains (California, USA) that exhibit dramatically different topographies despite their proximity and associated similar long-term climates. Informed by these measurements, we performed long-term and incremental alluvial fan experiments built by debris flows with systematically-varied fine particle content. We found that (1) decreasing fine particle content increases the variability of fan slopes and associated channelization dynamics, and (2) for all mixtures longer-term continuous alluvial fan experiments form more complex surface channelizations than repeated flows for the same total time, indicating the importance ofbothparticle sizes and timescales on alluvial fan surface morphology. 

This paper presents a discrete element model for simulations of the compaction process of hot mixed asphalt (HMA). The model is anchored by the concept of a fine aggregate matrix (FAM), which consists of the binder and fine aggregates. In the simulation, the coarse aggregates are explicitly modeled as composite particles. Meanwhile, the FAMis considered as the thick coating of the coarse aggregates with complex constitutive laws. Interparticle interactions include influences of (1) particle properties via Hertz–Mindlin relations; and (2) FAM properties via lubrication relationships. The lubrication relationships include a variable for viscosity for which we derive normal and tangential rate-dependent forms using rheology theory of dense granular-fluid systems, verified reasonable for our systems with the discrete element simulations and experiments with FAM. We assimilate these elements into gyratory compaction simulations of HMA of different aggregate size distributions. We compare these with experiments and find that this model is capable of capturing the measured effects of grain size distribution on the overall compaction behavior of HMA. We conclude by highlighting the advantages of this discrete element model for HMA compaction problems. 

Abstract Modeling transport, erosion, and deposition of nonuniform sediment over temporal intervals that are short compared to those characterizing channel bed aggradation and degradation remains an open problem due to the complex quantification of the sediment fluxes between the bed material load and the alluvial deposit. Parker, Paola, and Leclair in 2000 proposed a morphodynamic (PPL) framework to overcome this problem. This framework is used here to model the dispersal of a patch of gravel tracers in three different settings, a laboratory flume, a mountain creek, and a braided river. To simplify the problem, (a) the bed slope, bedload transport rate, and bed configuration are assumed to be constant in space and time (equilibrium), (b) sediment entrainment and deposition are modeled with a constant step length formulation, and (c) the PPL framework is implemented in a one‐dimensional (laterally averaged) model. Model validation against laboratory experiments suggests that, as the transport capacity of the flow increases, the maximum elevation‐specific density of sediment entrainment may migrate downward in the deposit. The comparison between model results and field data shows that the equilibrium solution can reasonably capture tracer dispersal. The equilibrium model can also reproduce subdiffusion and superdiffusion of a patch of tracers in the streamwise direction, depending on the magnitude of the short‐term bed level changes. Finally, the average tracer elevation in a cross‐section decreases in time because particles that are buried deep in the deposit are only rarely reentrained into bedload transport.